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Abstract 

We study CP-violating effects in neutralino production and subsequent decay 
within the Minimal Supersymmetric Standard Model with complex parameters M\ 
and \i. The observable we propose is a T-odd asymmetry based on a triple product 
in neutralino production e + e~ — ► XiTGi * = l>--->4, with subsequent leptonic 
three-body decay X2 ~* X?^ + ^~) £ = e,[i, a t an e + e~ linear collider with y^ = 
500 GeV and polarised beams. We provide compact analytical formulae for the cross 
section and the T-odd asymmetry taking into account the complete spin correlations 
between production and decay. We give numerical predictions for the cross section 
and the T-odd asymmetry. The asymmetry can go up to 10 %. 



1 Introduction 

In the Standard Model (SM) the phase in the Cabibbo-Kobayashi-Maskawa (CKM) ma- 
trix is the only source of CP violation. The small amount of CP violation in the SM, 
however, is not sufficient to explain the baryon-antibaryon asymmetry of the universe 
[1]. The Lagrangian of the Minimal Supersymmetric Standard Model (MSSM) contains 
several complex parameters, which can give rise to new CP- violating phenomena [2]. 
In the neutralino sector of the MSSM two complex parameters appear, which lead to 
CP-violating effects in neutralino production and decay. These parameters are the U(l) 
gaugino mass parameter M\ and the higgsino mass parameter \x. The phase of the SU(2) 
gaugino mass parameter Mi can be eliminated by an [/(l)-symmetry of the model. 

One of the main goals of a future e + e~ linear collider will be a careful study of the 
properties of supersymmetric (SUSY) particles [3]. The neutralinos x°, i — 1, • • • , 4 , will 
be particularly interesting, because the lightest neutralino Xi is expected to be the lightest 
SUSY particle (LSP), which is stable if i?-parity is conserved. The second lightest neu- 
tralino, X2, will presumably be among the lightest visible SUSY particles. Therefore, the 
study of production and decay of the neutralinos xl, [3, 4, 5] and a precise determination 
of the underlying supersymmetric parameters Mi, M 2 , \x and tan/5 including the phases 
0Mi and M of Mi and \x will play an important role at future linear colliders. The phases 
of M\ and \x cause CP-violating effects already at tree level. Therefore these effects could 
be large and thus be measurable at a high luminosity e + e~ linear collider. Methods to 
determine these parameters in neutralino and chargino production have been presented 
in [6, 7, 8, 9] . In particular the method of [9] is based on the analysis of the masses and 
production cross sections of only the light neutralinos Xi and X2 an d the light chargino xf 
and will allow the determination of the parameters \M\\, Mi, \\i\ and the phases 4>M t an d 
(pfj,. All these methods involve only CP-even quantities, in which the signs of the phases 
(f>Mi an d M can not be determined. For an unambiguous determination of <pM 1 and ^ 
including their signs one has to rely on CP-sensitive observables. 

The phases of the complex parameters are constrained or correlated by the experimen- 
tal upper limits on the electric dipole moments of electron, neutron and the atoms 199 Hg 
and 205 T1 [10]. In a constrained MSSM the restrictions on the phases can be rather se- 
vere. However, there may be cancellations between the contributions of different complex 
parameters, which allow larger values for the phases [11]. For example, in a constrained 
MSSM and if substantial cancellations are present, M is restricted to |0 M | < O.I71", whereas 
0Mi and (ft a, the phase of the trilinear scalar coupling parameter, turn out to be essentially 
unconstrained, but correlated with M [12]. Moreover, the restrictions are very model de- 
pendent. For example, when also lepton flavour violating terms are included, then the 
restriction on M may disappear [13]. Therefore it is necessary to determine the phases in 
an unambiguous way by measurements of CP-odd observables. 

A useful tool to study these CP-violating effects are T-odd observables, based on triple 
products of momenta or spin vectors of the particles involved [14, 15]. In this paper we 
study a T-odd asymmetry in neutralino production 

e + e-^xlx°2, k = l,---,4, (1) 



with subsequent leptonic three-body decay 

xl - xTr, (2) 

with £ — e, ii. Using the triple product of the initial electron momentum p e - and the two 
final lepton momenta pi+ and pg- , 

T = Pe+ ■ (pe- xp e -), (3) 

we define a T-odd asymmetry 

_ fsign{O r }|T| 2 rflips 
T /[Tulips ' [ } 

where f|T| 2 <ilips is proportional to the cross section cr(e + e~ — > XkX2 ~^ xtx®^^)- 
It is essential to include the spin correlations between production and decay, because 
otherwise the numerator and hence the asymmetry At would vanish. This observable 
has the advantage that it is not necessary to reconstruct the momentum of the decaying 
neutralino. The asymmetry At, Eq. (4), is odd under the naive time-reversal operation. 
By CPT it is a CP-sensitive asymmetry, if final-state interactions and finite-widths effects 
are unimportant. We will neglect these effects, because they are of higher order. By its 
definition, the asymmetry At is the difference of the number of events with the final 
lepton £ + above and below the plane spanned by pg- x p e -, normalised by the sum of 
these events. 

In this paper we will first present the complete analytic formulae at tree level for the 
cross section and for the asymmetry (4) of the processes (1) and (2), taking into account 
the full spin correlations between production and decay of the neutralino \2 [16]. We will 
study the dependence of the asymmetry and the cross section on the phases M and 0Mi 
and present detailed numerical results for an e + e~ linear collider with \/s = 500 GeV and 
polarised beams. Analogous CP asymmetries in neutralino production and subsequent 
two-body decays have been studied in [17]. CP asymmetries based on triple products 
in decays of scalar fermions have been discussed in [18]. CP-odd observables involving 
the polarisation of the outgoing r leptons from neutralinos decaying via two-body decays 
X2 -^ T-y t t — > Xi T±1 ~ T have been analysed in [19]. A systematic study of the impact 
of the supersymmetric phases on neutralino, chargino and selectron production at linear 
colliders has been performed in [12]. A Monte Carlo study of the asymmetry (4) for 
e + e~ — > X1X2 ~~ > XiXi^ + ^~) / = e, /x, including SM backgrounds and detector effects has 
been given in [20], however, no analytic formulae for the cross section and the asymmetries 
have been given there. In [21] the impact of the phases <pA T , 4>A t , <pA b , <fin an d $m x on 
the two-body decays of the third generation sfermions has been analysed in detail. The 
influence of these phases on the polarisation of the outgoing fermions in third generation 
sfermion decays has been studied in [22]. 

In Section 2 we shortly present the basics of the spin density matrix formalism. In 
section 3 we give the analytic expressions for the asymmetry At and the cross sections. 
Section 4 contains numerical results and discussions. In section 5 we present our con- 
clusions. Appendices A and B contain some details of the formalism and of the analytic 
expressions. 



2 Formalism 

2.1 Lagrangian and couplings 

The production process e + e~ — ► XiX2i * = 
(Fig. 1). In the decay process x° — * Xi 
(Fig. 2). 



, 4 , proceeds via Z° and e^ exchange 
e, // , Z° and £l^ exchanges contribute 



e+(p 2 ) 



e (pi) 




X°(P3) e+(p 2 ) 



eL,/? 



Xj(P4 





e (pi) X"(P4) e (pi) X 4 °(P3) 

Figure 1: Feynman diagrams of the production process e + e~ — ► x^X?- 




v~/0 _^ -zflp+t 



Figure 2: Feynman diagrams of the three-body decay Xi ~ ¥ xi 

The interaction Lagrangians for these processes are [23] (in our notation and conven- 
tions we follow [16]) 
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z°e+e- 






fj 
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2 cos 6vk 
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(5) 
(6) 



£«*2 = gfkiPRxlh + gfz R k £PLX°J R + h.c, (7) 

with m,n,k = 1,- ■ ■ ,4 and the couplings 

/£ = -v^— ^— (T 3 ^ - e, sin 2 9^)iV fe2 + e £ sin 9 w iV fel ] , (8) 

cos Wvk 

/* = -V2e e sme w [ttmG w NZ 2 - N* kl ] , (9) 

O'L = ~\{N m3 K 3 - N m ,N* n4 ) cos 2/3 - ^(N m3 N* A + N m4 N* n3 ) sin2/3, (10) 

qIIR _Q»L* , u \ 

^mn ^ van i \ ±± J 

L e = T 3e - e e sin 2 Q w , R e = -e £ sin 2 Q w , (12) 

where Pl,r = |(1 =F 7s), 5 is the weak coupling constant (g = e/cosGV, e > 0). e^ and 
T 3 £ are the charge (in units of e) and the third component of the weak isospin of the 
fermion £. Qyy is the weak mixing angle and tan/3 = v 2 /vi is the ratio of the vacuum 
expectation values of the Higgs fields. The unitary (4 x 4)-matrix N m k which diagonalises 
the complex symmetric neutralino mass matrix is given in the basis (7, Z, H®, H®) [4]. 

2.2 Amplitudes 



■~i\i Ao 



From these interaction Lagrangians the helicity amplitudes T p l 3 (a) and T D ^(a) for 
production and decay can be calculated. Here a denotes the production or decay channel 
and the respective exchanged particle. \ and Xj are the helicities of the neutralinos x® 
and x°j- I 11 TD,\A a ) we suppress the label of the helicity of the final neutralino xt- The 
amplitudes T p l J (a) for the production process e~(pi)e + (p 2 ) — »• X^ip^X^iPi) (see Fig. 1) 
are 

T A,A, = T A,A, ( ^ z) + T A,A 3 ( ^ gi) + T A,A, ( ^ ^ + ^ ^ ~j + ^ ^ g ^ ^ 

where 

T^(s, Z) = —^—A a (Z)v(p 2 )'f(L e P L + R e P R )u( Pl ) 

COS^ Ww 

«(p 4 , \ 3 )l»(0"±P L + 0'; t R P R )v(p 3 , X { ), (14) 

T^(i, e L ) = -g 2 f^f e L ;A t (e L )v(p 2 )P R v(p 3 , \)u[p A , A,)P L u( Pl ), (15) 

T? Xi (t, e R ) = -g 2 f^A\e R )v{p 2 )P L v{p 3 , Ai)«(p 4 , A,)P^( Pl ), (16) 

T^(«, e L ) = g 2 f^f^A u (e L )v( P2 )P R v( P ,, X 3 )u(p 3 , Ai)P L «(pi), (17) 

7^(u, e fl ) = g 2 f^f^ u {e R )v{p 2 )P L v{p^ X 3 )u(p 3 , X t )P R u( Pl ) (18) 

with 

A S (^) = 2 ' r , (19) 



and the Mandelstam variables 

s = (pi + p 2 ) 2 , t = (pi - Pi) 2 , u = (pi - p 3 ) 2 . (21) 

A S (Z), m z , T^, A t ' u (e i ^ R ), rrie LR , ^e LR denote the corresponding propagator, mass 
and width of the exchanged particle. The amplitudes T D> \.(a) for the decay x^ips) ~ ¥ 
X°k(P5)£ + (P(>)£~(P7) (see Fig. 2) read 

Tdx = T DX ( Si , Z) + T DX (U, l L ) + T DX (U, l R ) + T DX ( Ui , I L ) + T DX {u h l R ) (22) 

with 

T DX ( Si , Z) = 9 ——A s *(Z)u(p 7 )Y(LeP L + R e P R )v(p 6 ) 

cos^ fcV 

u{p,)l,{OlfP L + Ol*P R )u(p 3 , A,), (23) 

T DX {tJ L ) = -g 2 f e L jtA u ChMP7)PRv(p 5 )v(p3,\ l )PLv(p (i ), (24) 

W^*) = -9 2 fPkft ^{I R )u{p 7 )P L v{p,)v{p 3) \ t )P R v{p,), (25) 

T d ,a,(^,4) = +9 2 fkft A u '(IlM P7 )P r u( P3 , XMP^PMp,), (26) 

T DX {uJ R ) = +g 2 f«f«*A Ui (£ R )u(p 7 )P L u(p 3 ,\ i )u(p 5 )P R v(p 6 ) (27) 

and the kinematic variables 

Si = (pe + P7) 2 , *» = (P3 - Pe) 2 , Mi = (P3 - P7) 2 - (28) 

The amplitude for the whole process e + e~ — > x?Xj ~~ *■ Xi^ + ^~Xjj with x° decaying, can 
be written as 

T = A(x°)^T^T DA , (29) 

A* 

where A(%°) = l/[p§ — m 2 + otijFj], m^, Tj are the propagator, mass and width of the 
decaying neutralino xl- 

2.3 Cross section 

Following the formalism of [16, 24], the amplitude squared \T\ 2 of the combined processes 
of production e + e~ — > xlxl an d decay xl, ~* Xi^ + ^~ of the neutralino Xii with neutralino 
X°j being unobserved, can be written as 

m 2 = 4|A(x°)i 2 |p(x?x5)j>(x?) + E^x^oa j ' ( 3 °) 

where a = 1, 2, 3 refers to the polarisation state of the neutralino Xi, which is described by 
the polarisation vectors s a (x®) given in Appendix A. P{xlxT) an d D(x,i) are the terms of 
production and decay independent of the polarisation of the decaying neutralino, whereas 
T, P (Xi) an d Y/piyxl) are the terms containing the spin correlations between production and 

6 



decay. According to our choice of the polarisation vectors s a (x°), Eqs. (78) - (80), Ep/ P is 
the longitudinal polarisation, Ep/ P is the transverse polarisation in the production plane 
and Ep/P is the polarisation perpendicular to the production plane of the neutralino %°. 
The explicit expressions for these contributions are given in Appendix B. The complete 
expressions for the amplitude squared including spin-spin correlations for the case when 
both neutralinos are decaying and for longitudinally polarised e^ beams are given in [16]. 
The differential cross section in the laboratory system is 

I 7 

da = 8& l T ! 2 ( 27r ) 4( *V + P 2 - Z>i)<*Mp 3 • ' -^)> (31) 

b j_4 

where E^ is the beam energy and (flips (p^ ••• p?) is the Lorentz invariant phase space 
element. Integrating over all angles gives the cross section for the combined process of 
production and decay 

a = a (e+e- -> x°X?) • BR{% - xTH 

= iL /nx^ )^(x°)|A(x°)| 2 (27r) 4 5 4 (p 1+ p 2 -^^)dipsb3---P7). (32) 

b J i=4 

3 Triple product and CP asymmetry 

In this section we will derive the analytic formulae for the cross section a(e + e~ — ► XiXl ~~ *■ 
X?Xi^ + ^~) an d for the asymmetry At, Eq. (4). As can be seen from the numerator of 
At, for this purpose we have to identify those terms in \T\ 2 , Eq. (30), which contain a 
triple product of the form Eq. (3). Triple products follow from expressions ie^ upa k fl p u q p s' T 
in the terms Ep(x°) and EJ,(x°), where k, p, q, s are 4-momenta and spins of the par- 
ticles involved (see Appendix B, Eqs. (91) - (99) and (108) - (116)). The expressions 
it[ivpak tJ 'p u q p s a are imaginary and when multiplied by the imaginary parts of the respec- 
tive couplings they yield the terms which contribute to the numerator of At, Eq. (4). 
Hence only the second term of Eq. (30), which is due to the spin correlations, contributes 
to the numerator of the asymmetry At- It is convenient to split S P (x^) and SJ,(x°) into 
T-odd terms Ep (%°) and E^ (xl) containing the respective triple product, and T-even 
terms Ep (xD an d E^ (x°) without triple products: 

E P (x ? °) = E£°(x°) + E p ' E (x°), E a D (x°) = E£°(x?) + ^ E (x°). (33) 

Then the terms of |T| 2 , Eq. (30), which contribute to the numerator of A T are 

iti 2 d 4|A(x°)i 2 ^ [zfm^m + ^ a P E m^ a D °m] , m 

a=l 

where the first (second) term is sensitive to the CP phases in the production (decay) of 
the neutralino x°i- I n the following we derive explicitly the T-odd contributions to the 
spin density matrices of production and decay. 



3.1 T-odd terms of the production density matrix 

As can be seen from Eq. (99) in Appendix B, in Sp(x°) the expressions f% vanish for 
a = 1,3, because they do not contain three linearly independent vectors, hence E p ' (xl) — 
E p ' (xi) = 0- Only for a = 2 we get a T-odd contribution to the production density matrix 
Sp (x^), which is related to the polarisation of the neutralino Xi perpendicular to the 
production plane (see the expression for s 2 (x°), Eq. (79) in Appendix A). We obtain from 
Eqs. (92) - (95) for a = 2 



s p (Xi) 



AO, 



,2,0, 



^2,0 



2,0/ ~ ~ 



-,2,0/~ ~ 



Ep' (ZZ) + Ep' (Ze L ) + E p ' (Ze R ) + S P ' (e L e L ) + E P ' (e R e R 



(35) 



with 



Y?f{ZZ) 



-,2,0/ v ~ 



v^2,0/~ ~ > 

Ep (e L e L ; 



E^°(Ze R ) 



v^2,0/~ ~ \ 

Ep (epep) 



<f 



cos 4 W 



|A*(Z)| 2 (p 2 (l + Pe-)(l ~ Pe+) ~ L 2 e (l ~ P e -)(1 + Pe 



2(i?eO J7 L )(/mO J7 L )i/ 4 2 



-L £ (l-P e -)(1 + P e+ ) 



2 COS 2 ©14/ 

xPe{ A'(Z) [/^/^Oy*A**(e L ) - /£#C>yA"(g L )] / 2 } ; 
^(1 - P e -)(1 + P e+ )Re{A u (e L )A t *(e L )(ft) 2 (f" 3 ) 2 fl}, 

-P £ (l + P e -)(1-P e+ ) 



2 cos 2 614/ 

xPe{ A*(Z) [/* fljO^A^en) - f^O'^A u \e R )] / 4 2 }, 
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1 + P e -)(1 - P e+ )Pe{A«(gp)A'*(gp)(/ £ f ) 2 (^) 2 /|} ; 



(36) 

(37) 
(38) 

(39) 
(40) 



where P e - and P e + denotes the degree of longitudinal polarisation of the electron beam 
and positron beam, respectively. The function f%, Eq. (99), for a = 2 reads 



fl 



imjewpcPfiPiS^Ps. 



(41) 



Note that /| is purely imaginary. When inserted, for example, into the expression for 
Ep (ZeiS), Eq. (37), /| is multiplied by the factor i • Im{f^ i f^ J *O i j J }, which is non- 
vanishing only if the couplings are complex. Hence it gives a CP-sensitive contribution to 
the asymmetry At, which depends on the phases M and 0Mi- Analogous contributions 
come from the other terms in E p ' , Eqs. (36), (38) - (40). We have to multiply E p ' in 
Eq. (34) by E^ , for which we obtain from Eqs. (108) - (112), Appendix B, 



-,2,E/~0\ 



-i2,E, 



^2,E 



2,E, 



-<2,E, 



-<2,E, 



E^(xD = E^(ZZ) + ^(Z£ L ) + T^{Zi R ) + ^(hh) + E^W* 



(42) 



with 



Z 2 f(ZZ) = 8 



9' 



cos 4 Qw 



\A S \Z)\\R]-L] 



[(ReO^r-ilmOl^gl + lO'^igf-g!) 



^2,E 



VF(ze L ) 



^(Lh) 



X 2 l?(ZIr) 



V 



cos z <dw *■ L 

+ ftfe L k & Ui VL)(20'^g 2 2 + OM 

2^[|/^ri/^| 2 [|A^(^)| 2 ^ 2 - |A*<(I L )| 2 ^] 



V 



EBW*) = 2 5 



ifcite A-(Z) fX*^(£ R )(20 k fgi-O k r 9 i) 

COS z fc% *■ L 

+ ftfuA^R^ - X>** d - O'ifgl) 

\fZ\ 2 \f? k \ 2 \-\^{I R )\ 2 gl + |A'*(4)M 



where 



9l 
9l 
9l 



mi(p$P7)(P(>s 2 ), 

m{p<oPo){p 7 S 2 ), 
m k [(p3P 6 )(p 7 S 2 ) - (P3P7)(P6S 2 )}. 



(43) 



(44) 



(45) 



(46) 



(47) 



(48) 
(49) 
(50) 



The kinematic functions g 2 , g 2 , g 2 are real. When multiplied by the purely imaginary 
/|, Eq. (41), this leads to triple products sensitive to the CP phases in the production 
process, which in the laboratory system read: 



9l 


ft 


9l 


ft 


9l 


ft 



i2E b m i m j (p 5 p 7 )p 6 (p 1 x p 3 ), 

i2E b m i m j (p 5 p (i )p 7 (p 1 x p 3 ), 
i2E b m j m k {(p 3 p 6 )p 7 (p 1 x p 3 ) - (p 3 p 7 )p 6 (pi x p 3 )}. 



(51) 
(52) 
(53) 



As outlined above, these expressions will be multiplied in Eqs. (36) - (40) by the factors 
i ■ I'm{f[ i f[ :J *O ij L } etc. and contribute to the first term of Eq. (34) and, hence, to the 
numerator of the asymmetry Ax, Eq. (4). 



3.2 T-odd terms of the decay density matrix 

The second term in Eq. (34) is sensitive to CP violation in the neutralino decay process 
(2). The expressions g%, Eq. (116), in S^ (xi) contain three linearly independent vectors, 
however, E p ' (xl) (Eq- (61)) vanishes because s 2 (x°) is perpendicular to the scattering 
plane. Therefore, only the terms £p (xV) for a = 1,3 have to be taken into account. 
The CP-sensitive terms of the decay density matrix following from Eqs. (108) - (112) for 



a = 1, 3, are 



-iO,Oz ~0\ 



ia,0, 



^a,0 



a,0, 



^a,0, 



ia,0, 



Y^{Xi) = Z a ^(ZZ) + T^(Zt L ) + T^{Zl R ) + ^(i L i L ) + T^W R 



(54) 



with 



E a f{ZZ) 



Z a f(Z£ L ) 



VfihtL) 



9' 



COS 4 <$>w 

4g* 



A S \Z)\\L\ - R\) 2Re(O k t)Im(O k t)ig a 4 



L e Re\A s *(Z) -ffjf:O k ^A^(£ L ) 



^a,0 
J D 



,0/ 



(zi 



R 



COS 2 &w 

+ ftf^O'^A^Ch 
2g 4 Re{(ft) 2 (fL) 2 ^ U (h)A^Ch)9t}, 



9t 



cos 2 Q w 



R e Re\A^(Z) ~fgf£0'£*A*(i R ) 



+ feTft k O k lA^(e R ) 
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J^rWn) = 2g*Re{(f e ry(f? k yAHt R )A^(£ R )g:\ : 
where for a = 1, 3 we have 

g\ = imke^paS^v^ViVl- 



(55) 

(56) 

(57) 

(58) 
(59) 

(60) 



Now g%, a = 1,3, is purely imaginary. When inserted, for example, in Eq. (56) it is 
multiplied by the factor i ■ I'm{fj J i f^ k *O k f*} 1 which depends on the phases M and </>Mi 
and contributes to S^ . Analogous contributions follow from Eqs. (55), (57) - (59). The 
corresponding T-even terms of the production density matrix also entering in Eq. (34) 
are obtained from Eqs. (92) - (95), 

££ E (x°) = ^ E (ZZ) + Hf{Z~e L ) + Zf(Ze R ) + ££ E (e~ L e L ) + ^ E (e R e R ), (61) 



where 



^a,E 



Y?f{ZZ) 



T,f{Z~e L ) 



9 1 



|A S (Z)| 2 (i? 2 (l + Pe-)0- ~ Pe+) - L?(l - Pe-)(1 + P e+\ 



cos 4 Qw 

x I \d[f \\ft - ft) - [{ReO'^f - (ImO'lffMl 

Lt(l - P e -)(1 + P e+ ) 



(62) 



9' 



2 COS 2 Qy/ 

xRe{A°(Z) [2rfjL*0: 3 L *A"*(e L )f? - 2f£f%0'gA*'(e L )f2 

+ [fkf^O^A u *{e L ) + ftfiA^A^ieLM] }, (63) 



Sp E (e L e L ) 



9 -(l - P e -)(1 + Pe+) f|/|| 2 |^| 2 (|A u (e~ L )| 2 A a - |A*(e L )| 2 / 2 a ) 



+ Re{(ftr(ft j r^ u (eL)A t *(e L )f 3 a } 



(64) 



jp (Ze 



R 



-R e (l + P e -)(1-P e 

cR rR*/^i"R* \u* ( ~ 



2 COS 2 Qy/ 

xRe\A s (Z) 2jfJ%Ol K *A u *{e R )K + 2fffgO'f *'*(<-■„)% 
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^a,E/~ ~ \ 9 



1 



[f%f5*0"«^*{e R ) + f^fp^^{e R m] }, (65) 

cR\2\tR\2( \\u(~ \\2ra , ul/- \|2 /-a\ 



Sr(e^) = j(l + ^e-)(l-Pef)[|/i?l1/^(-|A«(e fl )|^ + |A'(g a )|^; 



M(/« )V^Y^ u {e R )^{e R )r,} , (66) 



where 



A° = rmbwjipis*), (67) 

/ 2 a = mi (p m )(p 2 s a ), (68) 

/3 = ^[(PiP3)(p 2 s a ) - ( P 2P3)(piS a )]- (69) 

The triple products sensitive to the CP phases in the decay read in the laboratory system: 



^ ft ■ 9l = irnim k (p2P4){ - E b p 5 (p 7 x p 6 ) - E 7 p 5 (p 6 x pi) 

i=l,3 . 

+ E 6 p 5 (p 7 x pi) + £ 5 Pi (pV x p 6 ) | , (70) 

X] ^ ' #4 = «"li"lfc(piP4)| - E b p 5 (p 7 x p 6 ) + E 7 p 5 (p 6 x pi) 

- ^6Ps(P7 X pi) - £ 5 Pl(P7 X p 6 )|, (71) 



a=l,3 



5^ /a • 92 = irn j m k ^[{p 2 pz) ~ (piP3)]E b p 5 (p 7 x p 6 ) 
0=1,3 

+ [(P2P3) + (P1P3)] 



E 7 p 5 (Pe x pi) - E 6 p 5 (p 7 x pi) - £ 5 Pi(P7 x p 6 ) j. (72) 



As already mentioned, these quantities will be multiplied in Eqs. (55) - (59) by the factor 
2 • Im{ftife£Oki*} e ^ c - an d contribute to the second term of Eq. (34). 

We take into account only the contributions to the T-odd asymmetry At which stem 
from the complex couplings and neglect the contributions from the widths of Z, £l and 
£ R in the propagators Eqs. (19), (20), because they are of higher order. 

4 Numerical Results 

In this section we analyse numerically the CP asymmetry At, Eq. (4), and the cross 
section for reactions (1) and (2), at an e + e~ linear collider with ^Js = 500 GeV and 
longitudinally polarised e 1 * 1 beams. In particular, we investigate the dependence on the 
phases of the complex parameters Mi = \Mi\e l ^ M i and fi = \^\e 1 ^^ . For our numerical 
analysis we choose two scenarios A and B, such that m^o < m^o+m z o and m^o < m^ to 
prevent two-body decays of x°. The SUSY parameters and the masses of x°, i = 1, . . . , 4, 
and £l, ^r are given in Table 1. The decay widths of the neutralinos have been computed 
with help of the program SPheno [25] . 
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Scenario 


A 


B 


|Mi| 


150 


90 


M 2 


300 


192.7 


Ia*I 


200 


352.4 


tan/3 


10 


10 


m i L 


267.6 


267.6 


m ln 


224.4 


224.4 


(<pM 1 , M ) 


(0.5tt,0) 


(0.57r,0.57r) 


(0.5tt,0) 


(0.57r,0.57r) 




135.4 


139.9 


88.5 


89.7 


A2 


182.2 


182.8 


176.2 


180.3 


A3 


213.3 


214.9 


359.4 


360.7 


m xl 


334.8 


331.6 


377.1 


373.6 



Table 1: Input parameters |Mi|, M 2 , |/x|, tan/3, m^ and m~ t and the resulting masses 
m^o, i = 1, . . . ,4, for (0Mi,0ft) = (0.5-7T, 0) and (0.57T, 0.57r). The parameters M 2 , |/i| and 
tan/3 in scenario B are chosen according to the scenario SPSla in [26]. All masses are 
given in GeV. 



In scenario A the neutralinos x\i X 2 > X3 are mainly i?-higgsino, higgsino-S-VT 3 , 
higgsino--B mixtures, respectively, where the B component of the X2 nas a relatively 
strong <f> Ml dependence. In Fig. 3 we show the total decay width r(x 2 ) an d the branching 
ratio BR(x2 -^ Xi^ + ^~), ^ = e or /j in scenario A as a function of (f>Mi for M = 0. As can 
be seen, the CP-even quantities r(x°) and BR(x2 ~* Xi^ + ^~) depend quite strongly on 
<f>Mi- T(x 2 ) varies between 34 keV for 4>Mi = 0.257T, 1.757T and 4.2 keV for </>Mi = tt, and 
BR(x2 ~ ¥ Xi^ + ^~) between 6.4 % for (J)m 1 = 0, 2n and 2.6 % for (pM 1 = tt. The parameters 
M 2 , |/x| and tan/3 in scenario B are chosen like in the scenario SPSla in [26], with xl, X 2 > 
x\ being mainly B-, W 3 -, higgsino-like, respectively. In this scenario T(x 2 ) varies between 
80 keV for <p Ml = 0, 2vr and 9.9 keV for <p Ml = vr, and B(x 2 -> X?^ + between 1.9% 
for 4>m, = 0.47T, 1.67T and 2.9% for (pM 1 = tt. In both scenarios we take for the slepton 
267.6 GeV, m^ = 224.4 GeV and for the squark masses m„ fl = 597.6 GeV, 
600.5 GeV, mg = 602.9 GeV for the first and second 
587 3 GeV, m L = 615.7 GeV (A b = 1000 GeV) for the bottom 



masses m 
generation, and m~ b 



599.0 GeV and m~ dR 



61 



squarks. In most of the examples we take values for the phases </> M and </>Mi which are in 
agreement with the constraints from the electron and neutron EDMs. In order to show 
the full phase dependence of the asymmetry and the cross section, we also give some 
examples where we vary </> M and <J)m 1 in the whole range, relaxing the constraints from the 
EDMs. 
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Figure 3: (a) Total decay width r^) and (b) branching ratio BR(x2 ~~* Xi^ + ^ ), ^ = e 
or /j,, in scenario A of Table 1 for 0^ = 0. 



4.1 



e + e" 



~o-o 

A1A2 



First we discuss neutralino production e + e 



+,-- 



X1X2 with subsequent leptonic three- 



body decay x\ ~^ Xi^ + ^~- I n this case one obtains CP- violating contributions from both 
production and decay. In Fig. 4 (a) we show the CP asymmetry At as a function of the 
phase <f)M 1 for scenario A defined in Table 1 for M = and two centre of mass energies 
t/s = 350 GeV and y/s = 500 GeV. The beam polarisations are fixed at P e - = —0.8 and 
P e + = +0.6. The CP asymmetry attains the largest values of about 11% (—11%) for 
\fs = 500 GeV at 0m x ~ 0.27T (4>m 1 ~ I.871"). For y/s = 350 GeV, i.e. closer to threshold 
of the production, At reaches values of ±13.5 % at 4>m 1 ~ 0.2-7T, 1.87T. Choosing M = 0.17T 
instead of M = leaves the asymmetry At almost unchanged. Note that At does not 
have its largest values at cj>m 1 = f (<Ami = ^f ), where the imaginary parts of the respective 
couplings are maximal, but at smaller (larger) phases. The reason is an interplay between 
the 0Mi dependence of the cross section, shown in Fig. 4 (b), and the 0Mi dependence 
of the numerator of At, Eq. (34). This numerator is a sum of products of CP-odd and 
CP-even factors, which essentially have a sin 0Mi and cos 0Mi behaviour, respectively (for 
M = 0). We plot in Fig. 4 (b) the cross section a = a(e + e~ — > X1X2) ' BR(x2 ~^ Xi^ + ^~) 
for the production and subsequent decay process e + e~ — ► x?x° ~~ * XiX?^ + ^ _ (Eq. (32)), 
summed over £ = e, //. Note that also the cross section has a rather strong 0m x dependence. 
At 0Mi = 7i" the production cross section <r(e + e~ —> X1X2) nas a maximum, whereas the 
branching ratio BR^x? ~^ Xi^ + ^~) has a minimum (see Fig. 3 (b)) resulting in the dip of 
the cross section for the combined process of production and decay at <pM 1 = ^ ■ 

In Fig. 5 we show the contour lines of the CP asymmetry At in the 0Mi~</v plane in 
scenario A (Table 1) for \fs = 500 GeV and two sets of beam polarisations P e - = —0.8, 
P e + = +0.6 and P e - = +0.8, P e + = —0.6. In the scenario considered the 0m x dependence 
is stronger than that on M . In the case of P e - = —0.8, P e + = +0.6 the largest asymmetries 
\A T \ ~ 10 % are reached for <f> Ml ss 0.27T, 1.87T, rather independent of 0^ (Fig. 5 (a)). For 
P e - = +0.8, P e + = —0.6 (Fig. 5 (b)) the maximal values of \At\ ~ 7.5 % are attained near 
4>Mi ~ 0.37r, <f)^ pa 0.47T and <f>Mi ~ 1-7tt, 0^ ~ 1.6-7T. The sign flip of At between Fig. 5 (a) 
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Figure 4: (a) CP asymmetry At (Eq. (4)) and (b) cross section a(e + e 
Xi^ + ^~), summed over £ = e,/i, in scenario A of Table 1 for 0^ 
p e+ = +0.6 and ^ = 500 GeV (solid), y/s = 350 GeV (dashed). 



xlxl>BR{xl 



0, P e - 



-0.8, 



and (b) and the weak M dependence follows from the structure of the Xi^L,R^ couplings 
and the facts that m^ < m^ and that the magnitude of the B component of X2 depends 
rather strongly on <pMi- For right-handed (left-handed) electrons the B component (B 
and W 3 components) of Xi an d X2 contribute, which leads to a more pronounced (pM 1 
dependence than that on c/y 

In Fig. 6 (a) - (d) we show the contour lines of the CP asymmetry A? and of the 
cross section for e + e~ — ► X1X2) X2 ~^ Xi^ + ^~ ' 1 summed over £ — e, /i, in the M 2 -|/i| plane 
for \M 1 \/M 2 = 5/3 tan 2 9 W , <j> Ml = 0.5tt and M = 0, tan/3 = 10, m ih = 267.6 GeV, 
mi = 224.4 GeV, at ^Js = 500 GeV and two sets of beam polarisations, P e - = —0.8, 
P e + = +0.6 and P e - = +0.8, P e + = —0.6. For both polarisation configurations At is 
largest in the region |/x| ~ 240 GeV and M-i > 300 GeV with maximal values At ~ 7.5 % 
for P e - = -0.8, P e + = +0.6 (Fig. 6 (a)) and A T w -10% for P e - = +0.8, P e + = -0.6 
(Fig. 6 (b)). For P e - = —0.8 (P e - = +0.8) the main contributions to the asymmetry At 
come from the Z-£l (Z-£r) interference terms, therefore large asymmetries are attained 
if both neutralinos have significant higgsino and gaugino components. In Fig. 6 (c) and 
(d) we show the cross section cr(e + e~ — > X1X2 "~ > XiXi^ + ^~), summed over £ = e,/i. 
For P e - = —0.8, P e + = +0.6 (Fig. 6 (c)) the cross section is larger than about 10 fb 
in the region with maximal A T (|/i| w 240 GeV, M 2 > 300 GeV). For |/i| > 400 GeV 
and M 2 ~ 125 GeV the cross section reaches values larger than 40 fb, the asymmetry, 
however, is rather small, At ~ —2.5 %. For P e - = +0.8, P e + = —0.6 (Fig. 6 (d)) the cross 
section has values between 1 and 10 fb in the region where At is maximal (|/i| ~ 240 GeV, 
M 2 > 300 GeV). 

In principle, also for the case of hadronic neutralino decays into heavy quarks, x 2 — > 
X\bb and x 2 ~~ > X? c c, it will be possible to measure the asymmetry At- We have calculated 
At in the M 2 -|/i| plane and found values of roughly the same order of magnitude as those 
shown in Fig. 6 for the leptonic decay. The cross sections are larger by roughly a factor 
5 - 10 in the main part of the parameter region of Fig. 6 and approximately of the same 
order of magnitude as in the leptonic case for \fj,\ > 400 GeV. The experimental errors, 
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Figure 5: Contours of the CP asymmetry At (Eq. (4)) in % for e + e — > X1X2 with 
subsequent decay x° — > Xi^ + ^ m scenario A of Table 1 for -y/s = 500 GeV and (a) 
p e _ = -0.8, P e + = +0.6 and (b) P e - = +0.8, P e + = -0.6. 



however, are expected to be larger than for the leptonic X2 decays, because the distinction 
of the b and b or c and c charges and the resolution of the jet momenta are worse than in 
the leptonic case [27]. 

We also studied the \Mi\ dependence of the asymmetry At and the cross section, 
relaxing the GUT relation for \Mi\ with all other parameters as in scenario B of Ta- 
ble 1. In Fig. 7 (a) we show the T-odd asymmetry At as a function of 0Mi for four 
values of \M\\ = 90 GeV, 95 GeV, 100 GeV and 105 GeV. Here maximal asymmetries 
\At\ ~ 4% are reached. For \M\\ = 105 GeV At shows nearly two complete oscillations, 
whereas for \M\\ = 90 GeV there is only one oscillation with extrema at 0Mi = 0.757T and 
1.257T. In both cases this behaviour is caused by cancellations between the contributions 
Sp (x°)S^ (xi) from production and Ep (x°)E^ (x°) from decay in Eq. (34), which 
have opposite sign. For \M\\ = 90 GeV the contributions from the decay are dominant, 
whereas for \M-\_\ = 105 GeV both contributions are of similar magnitude, resulting in two 
oscillations. The corresponding cross sections are shown in Fig. 7 (c). In Fig. 7 (b) we 



show A T as a function of 0„ for 



>Mx 



0.25vr and three values of \M t \ = 95 GeV, 100 GeV 



90 GeV because in this case m^o > m^a + m z o 



and 105 GeV. (We do not discuss \M\ 

for 0.57T < (p^ < 1.77T.) For \M\\ = 100 GeV At exhibits the strongest variation with 
(f)^ with maximum \At\ ~ 3% for M m 0.2ir and minimum \At\ ~ for M ps -k. The 
corresponding cross sections are shown in Fig. 7 (d). 

Furthermore we have analysed the dependence of the asymmetry A T on the masses 
of the right and left sleptons, assuming m^ and m^ are free parameters, for <pM 1 — 0.57T 



and 



with the other parameters as in scenario A of Table 1. For ^J~s = 500 GeV 



and beam polarisations P e - 



-0.8 and P e + 



+0.6 we find that the asymmetry A T 
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Figure 6: (a), (b) Contours of the CP asymmetry At (Eq. (4)) in % for e + e 



x\xl 



X? 



X? 



and (c), (d) contours of the corresponding cross 
), summed over £ — e, /J, , in fb, respectively, for 



with subsequent decay x 2 

section a(e + e" -> x°x°) • £.R(x° 

tan/5 = 10, m lh = 267.6 GeV, m lR = 224.4 GeV, |Mi|/M 2 = 5/3 tan 2 9 W , <f) Ml = 0.5vr and 

M = with v^i = 500 GeV and (a), (c) P e - = -0.8, P e + = +0.6 and (b), (d) P e - = +0.8, 

P e + = —0.6. The points mark scenarios A and B of Table 1. The dark shaded area marks 

the parameter space with rn^± < 103.5 GeV excluded by LEP. In the light shaded area the 

analysed three-body decay is strongly suppressed because m^o > mz + m^o or m^o > rn~ iR . 
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Figure 7: (a), (b) CP asymmetry At (Eq. (4)) and (c), (d) cross section <r(e + e — ► 
X1X2) ' BR{j(2 ~ ¥ Xi^ + ^~)j summed over £ = e, /x , in scenario B of Table 1 with (a), (c) 
M = 0, (b), (d) Ml = 0.25tt and \M±\ = 90 GeV (solid), \M X \ = 95 GeV (dashed), 
|Mi| = 100 GeV (dotted) and |Mi| = 105 GeV (dashdotted) for P e - = -0.8, P e + = +0.6 
and x/i = 500 GeV. 



is larger than 6 % for rug ~ 200 GeV and mi > 150 GeV. For m^ ps 200 GeV and 
m i ~ 300 GeV the asymmetry is smaller than —6 %. This is due to the fact, that for 
these beam polarisations the contributions from £l exchange to the CP asymmetry are 
dominant. Outside of these narrow regions around m^ m 200 GeV or m^ rs 200 GeV 
the asymmetry is very small. For opposite beam polarisations P e - = +0.8, P e + = —0.6 
the £r contributions to A T are dominant and the sign of the asymmetry changes: for 



rrin 



in 



200 GeV and m ih > 250 GeV it is A T < -11 % and for 



m,- 



300 < m~ ( < 850 GeV it is A T > 6 %. 



*L 



200 GeV and 



4.2 



e + e 



~o~o 

A2A.2 



In the reaction e + e~ — ► X2X2 there are only T-odd contributions from the decay, i.e. 
from the second term in Eq. (34). There is no T-odd contribution from the production, 
because the amplitudes Eqs. (14) - (18) are real and the first term in Eq. (34) vanishes. 
We will give numerical results for the case where one neutralino decays leptonically and 
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the other decays hadronically. If both neutralinos would decay into lepton pairs of the 
same flavour, it may be experimentally difficult to distinguish the two lepton pairs coming 
from different neutralinos. In Fig. 8 we show the contour lines in the M 2 -\fi\ plane of the 
CP asymmetry At and of the cross section with one neutralino decaying leptonically, 
X2 ~^ x\^ + ^~i summed over l = e, //, and the other neutralino decaying hadronically, 
X2 ~^ X1QQ1 summed over q — u,d, s, c, b, for y/s = 350 GeV and P e - = —0.8, P e + = 
+0.6. The other parameters are tan/5 = 10, m^ = 267.6 GeV, mi = 224.4 GeV, 
\Mi\/M 2 = 5/3 tan 2 9 W , 4>Mi = 0.57T and 0^ = 0. For this process A T is considerably 
smaller than for e + e~ — > xHxHi because for production of a pair of equal neutralinos only 
the decay amplitudes contribute to the asymmetry. The largest asymmetries are obtained 
f° r I A 4 1 ~ 400 GeV and M 2 ~ 150 GeV, where Xi an d X2 have gaugino character. The 
maximum values are At rs 1.8 %. The cross section, however, reaches values up to 10 fb 
in this region. The hadronic branching ratio J2 q = u dscb BRijxH ~^ XiQQ) i s between 50 % 
and 95 % in the main part of the parameter region of Fig. 8 and between 30 % and 50 % 
f° r I A 4 1 ~ 400 GeV. For y/s = 500 GeV the asymmetries are smaller with maximum of 
0.7% for \fi\ > 400 GeV and M 2 « 150 GeV. 

We want to remark that it may also be possible to distinguish the two lepton pairs, if 
one neutralino decays to an e + e~ pair, the other one to a fi + fi~ pair. Similarly it may be 
possible to distinguish the two quark pairs, if one neutralino decays to a bb pair and the 
other one to a cc pair [27] or pair of light quarks. However, in the case of both neutralinos 
decaying into lepton pairs the event rate is expected to be very low, whereas in the case 
of the decay into quark pairs of different flavours the experimental errors are expected to 
be much larger. 

4.3 e+e- -► x%X% 

In Fig. 9 we show the contour lines of the CP asymmetry At in the 0Mi~</v plane in 
scenario A (Table 1) for yt~s = 500 GeV and two sets of beam polarisations P e - = —0.8, 
P e + = +0.6 and P e - = +0.8, P e + = —0.6. As in the case of X1X2 production (Fig. 5) the 
<Pm± dependence is stronger than that on M in this scenario, since the B component of 
the decaying particle x\ depends stronger on (f>M 1 - The largest asymmetries of \At\ > 9 % 
are reached in the case of P e - = —0.8, P e + = +0.6 (Fig. 9 (a)) for 0^ ~ 0.257T, 1.757T 
and 0.6tt < <^ < 1.4tt. In the case of P e - = +0.8, P e + = -0.6 (Fig. 9 (b)) the sign 
of A T flips and the largest asymmetries \A T \ > 7% appear for <p Ml ss 0.257T, 1.75tc and 

0.87T < ^ < 1.27T. 

In Fig. 10 we show the contour lines of the CP asymmetry At and of the cross section 
for the process e + e _ — ► X3X2 with the subsequent decay X2 ~* Xi^ + ^~ \ summed over 
£ = e,/i, in the M 2 -\n\ plane for tan/3 = 10, m lh = 267.6 GeV, m~ lR = 224.4 GeV, 
\M 1 \/M 2 = 5/3 tan 2 6 W , <f> Ml = 0.57T and 0^ = for <Js = 500 GeV and two sets of 
beam polarisations P e - = —0.8, P e + = +0.6 and P e - = +0.8, P e + = —0.6. \At\ has two 
maxima in the regions \/j,\ w 240 GeV, M 2 w 300 GeV and |/x| w 350 GeV, M 2 w 140 GeV, 
respectively, for both polarisation configurations. The maximum values are At > 10% 
for P e - = -0.8, P e + = +0.6 (Fig. 10 (a)) and A T < -8% for P e - = -0.8, P e+ = +0.6 
(Fig. 10 (b)). The cross section reaches values larger than 10 fb (6 fb) in the region with 
maximal asymmetry around |/x| w 240 GeV, M 2 m 300 GeV for P e - = —0.8, P e + = +0.6, 
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Figure 8: Contours (a) of the CP asymmetry A? (Eq. (4)) in % for e + e~ 
subsequent decay x 2 ~~ *■ Xi^ + ^~ an d (b) of the corresponding cross sectk 
X^X 2 )'BR{x2 ~^ Xi^ + £~)'BR(x2 ~^ XiQ<j)i summed over £ — e, /j, and q — u,d, s, c, b, in fb, 
respectively, for tan/5 = 10, m~ lh = 267.6 GeV, m iR = 224.4 GeV, \M 1 \/M 2 = 5/3 tan 2 9 W , 
4>m 1 = 0.57T and M = with y^ = 350 GeV and P e - = —0.8, P e + = +0.6. The dark 
shaded area marks the parameter space with m-± < 103.5 GeV excluded by LEP. The 
light shaded area is kinematically not accessible. 



Fig. 10 (c) (P e - = +0.8, P e + = -0.6, Fig. 10 (d)). In the region around |//| « 350 GeV, 
Mo ~ 140 GeV the cross section is < 1 fb. 

We have also calculated the branching ratio BR(x® — > Xi^ + ^~)i Z — e or A 4 ? °f x!$- 
In parameter regions, where X3 decays via a three-body decay and where it is especially 
difficult to distinguish the two lepton pairs coming from different neutralinos, BR(x 

X? 



r,o 



) is always smaller than about 7%. 
Furthermore we have calculated the CP asymmetry A? for the the process e + e" 



X4X2 with the subsequent decay x 2 ~^ Xi^ + ^ 



e, \x. For the associated production of 



x\ and x\ t he accessible parameter space at a linear collider with ^fs = 500 GeV is rather 
small and is approximately constrained by M 2 + \fx\ < 500 GeV. For the same parameters 
as in Fig. 10 we find that the CP asymmetry has maximum values \At\ ~ 6 %. However, in 
parameter regions with \At\ > 2 % the cross section cr(e + e~ — > X4X2) ' BR{x\ 
summed over £ — e, /j, , is always smaller than 1 fb. 



xi 



5 Conclusions 

In this paper we have studied CP violation in neutralino production and subsequent three- 
body decay processes. As CP-sensitive observable we have chosen a T-odd asymmetry 
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Figure 9: Contours of the CP asymmetry At (Eq. (4)) in % for e + e — » X3X2 with 
subsequent decay X2 ~ > Xi^ + ^ m scenario A of Table 1 for -y/s = 500 GeV and (a) 
p e _ = -0.8, P e+ = +0.6 and (b) P e - = +0.8, P e + = -0.6. 



based on triple product correlations between the incoming and outgoing particles. A 
non-vanishing value for the asymmetry indicates directly CP-violating effects. We have 
provided compact analytical formulae for the cross section of the whole process as well as 
for the T-odd asymmetry for longitudinal beam polarisations. It is necessary to include 
the full spin correlations between production and decay. The asymmetry can be directly 
measured in the experiment without reconstruction of the momentum of the decaying 
neutralino or further final-state analyses. In our numerical results we have chosen rep- 
resentative scenarios and have shown that this T-odd asymmetry can reach values up to 
13 % in some parts of the parameter space of the unconstrained MSSM, even for small 
CP-violating phases as may be indicated by the EDM constraints. This observable will 
therefore be an important tool in the search for CP-violating effects in the neutralino 
sector and the determination of the phases of the complex SUSY parameters. 
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Figure 10: (a), (b) Contours of the CP asymmetry Ax (Eq. (4)) in % for e + e~ — *■ X3X2 
with subsequent decay X2 ~^ Xi^ + ^~ an d ( c )> (d) contours of the corresponding cross 
section o"(e + e~ — >• X3X2) '-^-^(X2 ~~ ^ Xi^ + ^")j summed over £ = e, \i , in fb, respectively, for 
tan/3 = 10, m £ - L = 267.6 GeV, m~ iR = 224.4 GeV, \M 1 \/M 2 = 5/3 tan 2 9 W , <j> Ml = 0.5tt and 
M = with v^i = 500 GeV and (a), (c) P e - = -0.8, P e + = +0.6 and (b), (d) P e - = +0.8, 



P, 



—0.6. The dark shaded area marks the parameter space with m~± < 103.5 GeV 



excluded by LEP. The light shaded area is kinematically not accessible or in this area the 
analysed three-body decay is strongly suppressed because m^o > m z + m^o or m^a > mi , 
respectively. 
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A Kinematics 

We choose a coordinate frame in the laboratory system, where the momenta are 

Pi = E b (l, -sine, 0, cos6), (73) 

p 2 = E b (l, sine, 0, -cose), (74) 

p 3 = (^,0,0,-g), (75) 

p 4 = (E v 0,0,q). (76) 

6 is the scattering angle between the incoming e~(pi) beam and the outgoing neutralino 
X^iPi), the azimuth $ can be chosen equal to zero. The energy and the momenta of the 
neutralinos x'ji.Pi) an d X°i(Pz) are 

s + mf — m 2 s + m 2 — m 2 \/K s i m n' m j) 

Ei = -= l , Ej = 3 -= -, q = -¥ F , (77) 

2^ ' 3 2^7s ' H 2^1 K J 

where rrii, rrij are the masses of the neutralinos and A is the kinematical triangle function 
X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz. 

The polarisation vectors s^(x°) and s^(x°) (o, b = 1,2,3) of the neutralinos in the 
laboratory system are 

^(xh) = (0,Tl,0,0), (78) 

s^ixh) = (0,0,1,0), (79) 

s^ixh) = —(q,0ATEij), (80) 

m i,j 

where s 3 (Xi,j) describes the longitudinal polarisation, s 1 (x° J ) the transverse polarisation 
in the scattering plane, and s 2 (x° ,) the transverse polarisation perpendicular to the scat- 
tering plane. 

B Explicit expressions of the amplitude squared 

B.l Production 

Here we give the analytic expressions for the production spin density matrix. The terms 
P(x°i)(S) = P{ZZ) + P(Ze L ) + P(Ze R ) + P{e L ~e L ) + P(e R e R ) (81) 

are independent of the polarisation of the neutralinos [16, 28], 

P{ZZ) = -4?^\^(Z)\ 2 (Rlc R + L 2 c L ) 

COS 4 Ww 



x {|Oy| 2 (A + f 2 ) - [{ReO'^f - {ImO'g?\h 



P(Ze L ) = ^-^—L e c L 

2 cos 2 k3 w 



iS2) 

4 
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xRe{ A S (Z) [ - (A^^ftftAt + A u *{e L )f^O'^)h 
+ 2A t *(e L )ftrt j O'; j L f 1 + 2A u *(e L )f e L i ffrO'; j L *f 2 ] }, 

P(e L e L ) = ^{l^H^rdA^lVi + IA^)! 2 ^) 

-^{(/^) 2 (^) 2 A u (g L )A'*(g L )/ 3 }}. 
P(Ze R ), P{e R e R ) are obtained by the following substitutions in Eqs. (83), (84): 
L e ^R £ , c L ^c R , A^(e L ) - A*-«(e fl ), O'g^O^*, /£,-->/&. 
The beam polarisation weighting factors are given by 

c L = (1-P e -)(1 + P e+ ), 
c H = (1 + P e -)(1-P e+ ) 

and the kinematical factors in invariant form are 
h = (pm)(p2Ps), 

/ 3 = mimjipip 2 ). 
The terms 

S« (x°) = Z a P (ZZ) + £« (Ze L ) + Z% (Ze R ) + ££(e L e L ) + Z a P (e R e R ), 



(83) 

(84) 

(85) 

(86) 
(87) 



(89) 
(90) 

(91) 



which depend on the polarisation of the neutralino x? an d which are relevant for the 
calculation of the CP asymmetries, are given by 



^ a P (ZZ) 



Y, a p{Zei 



9 



COS 4 6 w/ 



\A s {Z)\\R}c R - L\c L ){\0"J\\r 2 - n 



9' 



[{ReO'^f - {ImO'^f\n + 2(ReO'g)(ImO" i f-)if2 



-L e c L 



(92) 



2 COS 2 Qy/ 

xRe{A s {Z)\2f^*O^A u *{e L )n - 2fi?%O^A»(e L )f2 
+ [fU%0'^A u *{e L ) + fg%0'lj"&{eL)\n 

+ \JiC fip"^ A'\e L ) - ftf^O'^A u \e L )]ft] }, (93) 



Z? P (e L e L ) 



9 A 



l^| 2 |^| 2 (|A M (e~ L )| 2 A a -|A'(e L )| 2 / 2 a ) 



+ Pe{(/^) 2 (^.) 2 A u (e i )A**(e L )(/ 3 a + / 4 a )} 
Yj°p{Ze R \ Y* P (e R e R ) are obtained by the substitutions Eq.(85) and 

/ 1,2,3 ~~ *■ ~~ J 1,2,3 an< ^ J A ~ ¥ J A 



(94) 



(95) 
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in Eqs. (93), (94), where the kinematical factors in invariant form are given by 

]l = rm(P2Pi)(pis a ), (96) 

ft = m t ( Pm )(p 2 s a ), (97) 

/ 3 a = m j [ipm){p2S a ) - ( P 2P 3 )(pis a )}, (98) 

ft = im^ vpa fy\s ap V l. (99) 

The polarisation vectors s a (x®) of the neutralino x° are given in the laboratory system 
by Eqs. (78) - (80). 

B.2 Leptonic 3-body decay 

Here we give the analytical expressions for the different contributions to the decay density 
matrix for the three-body decay x^ipz) ~^ X°k{.P5)^ + {.Po)^ {.Pi) ^ where we sum over the 
spins of the final-state particles [16, 28]. The contributions independent of the polarisation 
of the neutralino Xi are 



D(Xi) = D(ZZ) + D(Z£ L ) + D(Z£ R ) + D(£ L £ L ) + D(£ R £ R ) 



(100) 



with 



D(ZZ) = 8 



cos 4 Qw 



D(Z£ L 



D(£ L £ L 



0" k h\9i + 92) + [{ReO'^f - {ImO" k ff\g, , 

f e L Ji L k *^*Ch)(20^ gi + 0'^*g 3 ) 

+ ftf^ u '*(h)(20'^g 2 + O'ggs) 



A^-—L e Re{^(Z) 
cos 2 <d w I 



2g 4 



\fk\Vk\ 2 (\^ u Ch)\ 2 gi + \^Ch)\ 2 92) 

+ Re{(ft) 2 (fk) 2 ^'Ch)A^(£L)}9s 



(101) 



(102) 



(103) 



The quantities D(Z£ R ), D(£ R £ R ) can be derived from Eqs. (102), (103) by the substitutions 



L e ^R e , A u ' Ui (e L ) -> A u ' Ui (£ 
The kinematical factors are 



R 



o" L 



q"R 

w ki i 



f L — >■ f H 
JtLk JtLk- 



(104) 



9\ = (P5P7)(P3Pe), (105) 

92 = (P5P6)(P3Pr), (106) 

g 3 = mim k (p & p 7 ). (107) 

The contributions which depend on the polarisation of the decaying neutralino Xi are 



*£(*") = ^d(ZZ) + Z a D (Zh) + Vh(Ztn) + ^dWl) + TT D {t^ R ) 



108) 
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with 1 

^d(ZZ) = 8^— \A«(Z)\*(Rl-L*) 
cos 4 k3 w 

x " - [{ReO'^f - (ImO'£)*\£ + \0% | 2 (tf - g a 2 ) 

-2Re(0'£)Im(0'£)igi], (109) 

Vl>(Zh) = -%^URe{^{Z)\fUk^ U \h){-20^g a l +0'^{gt-gt)) 
cos^ fc% *• L 

+ fufk^ih^O'^g* + Cff (<£ + ^))] }, (110) 
X£&I L ) = 2( 7 4 [|/^| 2 |/£| 2 [|A , "(^)r^-|A*'(I L )|V] 



+ MU^V&Y&iiL^tiLM + s 4 a )} 



in: 



The contributions T l c j D (Z£ R ),T> c ]r ) (£ R £ R ) are obtained from Eqs. (110), (111) by applying 
the substitutions in Eq. (104) and in addition 

9i,2,s^-9i,2,3, 9l^9l (112) 

The kinematical factors are 

9l = m i (p 5 p 7 )(p 6 s a ), (113) 

g a 2 = m l (pr )P6 )(p 7 s a ), (114) 

gl = m k [{p 3 p 6 ){p 7 s a ) - {p3P 7 ){P6S a )}, (115) 

g\ = imkC^s^p^p*. (116) 
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